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Abstract 



We discuss the exotic hadron structure and hadron-hadron interactions in view of heavy ion collisions. First, we demonstrate that 
a hadronic molecule with a large spatial size would be produced more abundantly in the coalescence model compared with the 
statistical model result. Secondly, we constrain the AA interaction by using the recently measured AA correlation data. We find 
■ that the RHIC-STAR data favor the AA scattering parameters in the range 1 /a < -0.8 far 1 and r eff > 3 fm. 

© 2012 Published by Elsevier Ltd. 
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1. Introduction 

The first dozen of years in the 2 1 st century may be recognized as the starting point of the exotic hadron renaissance. 
In textbooks, hadrons are explained as qq (mesons) or qqq (baryons) composites, and many of the hadron masses are 
well described in the quark model. This common understanding of hadrons becomes doubtful in these years. Starting 
from D J y(2317) Jl]], we have found many hadronic states which we cannot understand in the naive quark model. For 
example, a penta quark state ududs is claimed to be observed at LEPS |0], while its existence is still controversial fl. 
Z + (4430) is a typical and clear exotic hadron [4]: Its mass is close to D\D* threshold and we expect it contains cc, 
and it has a positive charge. Thus the minimum quark content of Z + (4430) is ccclu. 

Understanding the structure of exotic hadrons is important in order to construct a new scheme over the quark 
model to categorize hadrons including normal and exotic ones. There are mainly two-types of structure considered 
for exotic hadrons. One of them is compact multi-quark structure, and the other is hadronic molecule structure. These 
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two types of structure would have different sizes. We expect that a multi-quark state has a similar size to normal 
hadrons, while the deuteron, a well-known hadronic molecule, has a much larger size than normal hadrons. 

Another aspect of exotic hadron physics is that it is related to the hadron-hadron interaction. The existence of 
A(1405) below the KN threshold leads to various aspects of KN interaction. Similarly, once the pole position of the 
S = —2 dibaryon (//) is fixed, AA interaction is strongly constrained. 

High energy heavy ion collisions would provide unique information on exotic hadron structure and hadron-hadron 
interactions. Various hadrons are produced abundantly in heavy ion collisions, and it is natural to expect that exotic 
hadrons are also produced. The dynamics of high energy heavy ion collisions is so complex that statistical argument 
becomes valid; We may dare to say, heavy-ion collisions are simple and clean. Via the hadron-hadron correlation 
measurement, it is in principle possible to extract the resonance pole above the strong decay threshold. 

In this proceedings, we discuss the exotic hadron structure and hadron-hadron interaction in view of heavy-ion 
collisions. In Sec. [2] we demonstrate that the production yield is sensitive to hadron size. In Sec. [3] we discuss AA 
correlation and its relation to AA interaction. 



2. Exotic Hadron Yields in Heavy Ion Collisions 

Several mechanisms have been proposed so far to gain energy in exotic states. Based on the diquark picture and 
color-magnetic interaction, exotic hadron states including heavy-quarks should exist in some channels. The color- 
magnetic interaction is proportional to l/mj/mj, where m, is the quark mass. Then in a state made of QQud (Q 
denotes a heavy quark), a diquark (ud)(QQ) component would be favored rather than a mesonic molecule component 
(Q u )(Qd)- Tl c ( J n = = 0) is one of theoretically proposed hadronic states made of ccud fl5[] . Since the strong 
decay to DD is forbidden by the angular momentum conservation and its predicted mass is lower than the D\D, T* c 
may have a small width. A hadronic molecular state with the same quantum number is also predicted in the pion- 
exchange model |6]. The pion couples D(0) and D*(l~) states and pseudo-scalar and vector meson masses are close 
with heavy-quarks (heavy-quark symmetry), then the two "states" (D*D and DD*) couple strongly and T* c can gain 
energy. 

It is a challenge to clarify the structure and mechanism for each of the exotic states to exist. One of the key 
quantity to distinguish a multi-quark state and a hadronic molecule state is the size: When it is a multi-quark state 
such as a diquark pair (ud)(QQ), the confining force between the diquarks would make the system compact. When it 
is a hadronic molecule such as a D*D bound state, the exotic hadron size would be determined by its binding energy 
or the range of the pion-exchange. Recent work on A(1405) sheds light on this idea |7]. They have proposed that the 
evidence of the KN picture of A(1405) may be found in the negatively large squared charge radius, which may be 
observed in the electric form factor. 

We have proposed that we can utilize high-energy heavy-ion collisions to obtain knowledges on the exotic hadron 
size J3l. At RHIC and LHC, abundant hadrons are produced and their yield ratio is well described by the statistical 
model, which assumes thermal equilibrium at freeze-out H . We expect the statistical model also works for the exotic 
hadron production, which is calculated to be frequent enough. One of the problems of the statistical model is its 
prediction power for resonance states. For example, the statistical model overestimates the yield of A(1520) by a 
factor of two or more. This discrepancy is explained based on the coalescence (recombination) model, which uses the 
internal wave function thus includes the angular momentum effects fioll . Since the reaction time at RHIC and LHC 
is not very long, it may not be reasonable to understand the equilibrium hadron production literally and we may need 
underlying hadron production mechanisms which result in statistical distribution of ground state hadrons. 

In order to discuss the production mechanism dependence, we compare the results of the coalescence model and 



the statistical model at RHIC and LHC. We adopt the setup proposed by Chen and his collaborators 111 111 . Hadron 
yields in statistical model, 



= ghVH r°° p 2 dp 

h ~ 2n 2 Jo y- l e ^-^l T " ± 1 



N, = £ £ (1) 



are obtained at the transition (hadronization) temperature T — Th — 175 MeV and volume V = Vh- gh, Jh and 
are the degeneracy, fugacity, and chemical potential of the hadron h. Light quarks (u, d and s) are considered to 
reach approximate chemical equilibrium. Since charm and bottom quarks are produced from initial hard scattering 
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Figure 1. Coalescence / Statistical model ratio i?^ s as a function of mass (left) and Tua (right). 



and their numbers are much larger than the equilibrium values, we put the fugacity y/, > 1 for hadrons which contains 
charm and bottom quarks. We have fixed the fugacity of hadrons containing charm and bottom quarks to reproduce 
the expected number of charm and bottom quark number from initial hard scattering. 
The yield of a hadron h in the coalescence model is given as, 



AT' =Sh 
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where / t h and fw are the Wigner function of the constituents and the intrinsic states, respectively, T and V are the 
temperature and volume at coalescence, gj, Nj and ntj are the degeneracy, yield and mass of the j-th constituent 
(hadron or quark), M is the sum of constituent masses, and l, denotes the orbital angular momentum of the /-th 
Jacobi coordinate and their sum is L. Coalescence model calculations are performed at the hadronization temperature 
T - Th and the freeze-out temperature T = Tp for the quark and hadronic coalescence, respectively. We have 
assumed a harmonic oscillator wave functions, whose frequency is given asw = 550 MeV for hadrons made of u and 
d quarks. For hadrons containing strange, charm and bottom quarks, we fit the statistical model results of A(1115), 
A c (2286) and A fc (5620) and obtain lo s = 519 MeV, u> c = 385 MeV and oo b = 338 MeV, respectively. 



In Fig. \T\ we show the ratio of the hadron yields in the coalescence and statistical models, R 



where jv™ dl / stat denotes the hadron yields per unit rapidity in the coalescence or statistical model. First, we note that 
the ratio for normal hadrons (open triangles) is in the range of 0.2 < R^ s < 2 (gray band). Here "normal" hadrons 
are defined as particle states considered to be made of qq and qqq for mesons and baryons, respectively; lowest mass 
states for given quantum numbers ( J 71 and flavor) of pseudoscalar mesons ( J n = 0~), vector mesons (1~), and l/2 + 
and 3/2 + baryons. We also categorize JV(1440)(l/2 + ), JV(1520)(3/2-), AT(1535)(l/2-) and £>i(2420)(l + ) as normal 
hadrons. Secondly, the coalescence model is found to predict smaller yields of compact multi-quark states. When we 
use the same hadron size parameter as that of normal hadrons, an addition of a s-wave, p-wave, or af-quark leads to a 
suppression factor of 0.36, 0.093, or 0.029, respectively fl. Thus compact multi-quark states are suppressed lfl2ll . 

Another interesting feature found in the coalescence-statistical ratios is the enhancement of spatially extended 
hadronic molecules. Let us consider the two-body .s-wave coalescence in isotropic environment. The coalescence 
yield is given as the convolution of the intrinsic Wigner function and the thermal distribution of the relative coordinate, 
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where Ax (Ap) is the width in the intrinsic Wigner function in the spatial (momentum) coordinate, p is the reduced 
mass, and T is the temperature at coalescence. We have assumed here the spatial Gaussian source with the radius 
R. When the minimum uncertainty AxAp = h/2 is assumed, the above yield shows a maximum when the spatial- 
to-momentum width ratio of the intrinsic Wigner function is the same as that of the source, which reads fuo - 
yJh 2 T/2fxR 2 , As an example, in the case of T - 170 MeV, \i = 500 MeV, R — 5 fm, the optimal value of the oscillator 
frequency is ha> = 16 MeV, which is much smaller than that of normal hadrons, hat = (300 - 600) MeV. In the 
right panel of Fig. Q] we show the hco dependence of R^ & . In the present calculation, we have assumed that the source 
size is large enough. These results include three-body hadronic molecules, but the trend is the same. The coalescence 
favors hadrons whose shape in the phase space is similar to that of the source, then the large source size and moderate 
T prefer extended hadrons in coalescence. 

If the coalescence is the underlying mechanism of the statistical model, the coalescence model would give better 
predictions of the hadron yields including resonances and exotic hadrons, and we can utilize high-energy heavy-ion 
collisions as a ruler of the hadron size; smaller and larger yields for compact multi-quark states and spatially extended 
hadronic molecule states, respectively. 



3. Exotic Interaction from heavy ion collisions — AA interaction — 

Where is the S = — 2 dibaryon, H ? This is a long standing problem in hadron physics. In 1977, Jaffe pointed 
out that double strange dibaryon made of 6 quarks (uuddss) may be deeply bound below the AA threshold due to 
the strong attraction from color magnetic interaction |fl3l . Dedicated experiments have been performed to find the H 
particle in these 35 years. Deeply bound H was denied by the observation of double A hypemuclei. For example, a 
double A hypernucleus AA He was found to decay weakly in the Nagara event, and the observed energy of AA He is 
6.91 MeV(= fi AA ) below the 4 He + AA threshold Q. If the mass of H is below 2M A - B AA , AA He should decay 
to 4 He + H strongly. The reason why the attraction is weaker than expected may be the determinant-type 3-quark 
interaction [15], which is repulsive in the H channel [16]. While the deeply bound H is denied, the attraction in the 
H channel may generate a pole in the weakly bound or resonance region. There are some hints in recent experimental 



and theoretical studies. The KEK-E522 experiment observed a bump in the AA invariant mass spectrum 01711 . Recent 
lattice QCD studies imply that H should exist as a bound state in the SU(3) limit and/or with heavy pion masses [18]. 
Thus the physics of the H particle is a long-standing as well as current problem. 

Existence of the H particle state is closely related to the AA interaction. The H particle pole is, if exists, near the 
AA threshold, and it is natural to expect that H contains a significant component of AA. AA interaction is important 
also for the dense matter equation of state (EOS). In many of theoretical calculations, A fraction in dense neutron star 
matter is compatible with the neutron fraction, then the strength of the AA interaction may affect the EOS. Until now, 
available information on AA interaction is scarce. We know that it is weakly attractive from the AA bond energy in 
AA He, AB AA = B AA ( AA He) - 2fi A ( A He) 0.6 MeV. From Aff AA ( AA He), the scattering length and the effective range 
in the AA L S channel are obtained as (a , r eff ) = (-0.77 fm, 6.59 fm) S or (a , r eft ) = (-0.575 fm, 6.45 fm) TO , 
but in principle we cannot determine two low energy scattering parameters from one observed value of AB AA . 

Thus other observational information on AA interaction has been desired. One of the ways is to observe the 
binding energies of various double A hypemuclei, as planned in the J-PARC E07 experiment. Another available 
observable is the AA correlation in nuclear reactions. Actually, KEK-E522 experiment fll7ll has demonstrated that 
AA invariant mass spectrum is enhanced in the low energy region compared with the phase space estimate and the 
classical transport model calculation, implying that AA interaction is attractive. In high-energy heavy-ion collisions, 
abundant A particles are produced, and we can measure the AA relative momentum correlation, which contains 
information on AA interaction. This idea is not new. It was proposed in '80s that we can fix resonance parameters, 
when the source size is small [21]. The correlation at low relative momenta was proposed to be useful to discriminate 
the sign of the scattering length ao, provided that the source size is large [22]; When AA has a bound state (ao > 0), 
the scattering wave function must have a node at r - ao in order to be orthogonal to the bound state wave function, 
then we may find the suppression of the correlation. Now RHIC and LHC have the vertex detectors, and we can really 
obtain the AA correlation data in heavy-ion collisions. 

We here discuss AA correlation in heavy ion collisions in view of AA interaction. Two particle correlation at low 
relative momentum from a chaotic source is known to be sensitive to the source size and the two particle interac- 
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tion 12311 . The AA correlation function is given as, 

J dx i dx 2 S (x , , p + q)S (x 2 , p - q) | ^ H (x 1 2 , q) 



Caa(<?) = 



J dxidx 2 S (xi,p + q)S(x 2 ,p - q) 
1 - I exp(-4^ 2 ) + I J <irS 12 (r) (^„(r)| 2 - \j (qr)\ 2 ) , (4) 



where if/^Hx, q) is the relative wave function having the relative momentum q in the final state, xo is the relative wave 
function in the i-wave, and S denotes the source function. In obtaining the second line, we have made following 
two approximations. (1) The single particle source function has a Gaussian profile whose width is independent of 
momentum, then the source function in the relative coordinate is given as S i 2 (r) = (2R -00 3 exp(-r 2 /4/? 2 ). (2) Only 
the s-wave relative wave function is modified from the free case by the AA interaction. 
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Figure 2. Scattering parameters from AA interaction models. 



The correlation function is determined by the source function and the relative wave function, the latter of which 
is sensitive to the AA interaction. We compare the results with several types of AA interactions. The first type of AA 
interactions is the Nijmegen models Il24ll . which are based on a meson and meson-pair exchange picture of baryon- 
baryon interactions. The second type of baryon-baryon interaction is the quark model interaction fss2 N25I1 . which 
takes account of the Pauli blocking at the quark level, gluon exchanges between quarks, and meson exchanges. We 
also compare the results of a one-boson exchange AA interaction, Ehime potential [26]. Since the Ehime potential is 
proposed before the Nagara event, it assumes a smaller AA bond energy, AZ?AA( AA He) = 3.6 MeV, than that obtained 
in the Nagara event. In actual calculations, we use two range Gaussian potentials which fit the scattering length and 
effective range for Nijmegen and Ehime potentials. For fss2, we use a phase-shift equivalent local potential, derived 
by using the inversion method based on supersymmetric quantum mechanics 125l 12711 . Scattering parameters (ao, r e ff) 
of these interactions are shown in Fig. [2] 

Now we shall try to extract the information on the AA interaction as well as the source size from the AA correlation 
data. In the left panel of Fig. [3] we show the source size dependence of the AA correlation C(q). Calculated results 
using fss2, as an example, are compared with the RHIC-STAR data lEijll . When the source size is small, C(q) - 1 has a 
long tail in the high momentum region. Thus we can obtain the source size by fitting the correlation at high momenta. 
The optimal source radius R depends on the AA interaction, but the data suggest that the apparent source size of A 
is around R ~ 2 fm, which is smaller than the pion and kaon source. This difference may be due to the flow effects. 
In the right panel of Fig. [3] we show the size dependence of the AA correlation with flow effects. Transverse flow 
generally extends the range of finite correlation to a higher momentum region. Extending the correlation to higher 
momentum means that the flow makes the apparent radius smaller than the actual source size. We find that the A 
source size may be in the range 3 fm < R < 4 fm when we take a reasonable flow parameter rjf — 0.5 (the transverse 
rapidity is given as Yj = rifrj/R)- This source size would be consistent with the proton source size. 
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Figure 3. Left: Source size dependence of the AA correlation. Right: Calculated and observed AA correlation. 



After fitting the tail region, we can discriminate the AA interaction from the behavior of C(q) at small q. In the 
left panel of Fig.|H we compare the results from several AA interactions. It seems that Nijmegen model D (ND) with 
the hard core radius of R c = 0.56 fm, Nijmegen soft-core 97 model (NSC97f), and quark model interaction (fss2) are 
consistent with the RHIC-STAR data. It should be noted that these results are obtained in a simple setup; the single 
channel calculation, no feeddown effects, and no flow effects. 
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Figure 4. AA correlation obtained by using several AA interactions. The source size is chosen to fit the high momentum tail region of the 
correlation. Left (right) panel shows the results without (with) I. decay effects. 



For quantitative discussions, we need to consider the feeddown effects from heavier particles. The feeddown 
effects have been known to be decisive for the pp correlation. The pp correlation at low momentum is suppressed 
by the Coulomb repulsion, and the Gamow factor correction recovers C(q) to be around unity. In high-energy heavy- 
ion collisions, we have protons also from the decay of heavier particles, such as A — > pn~. There is no Coulomb 
suppression in the pA channel, and the pionic weak decay does not change the baryon momentum much. As a 
result, pA correlation in the reaction region strongly affects the pp correlation. Compared with the pp correlation, the 
feeddown effects on the AA correlation are expected to be small. There is no Coulomb suppression in the AA channel, 
and the particles which decays into A are limited. Furthermore, it is possible to exclude the A from weak decay such 
as H~ — > An~ using the vertex detectors, if necessary. There exists an exception, E° — > Ay, which we cannot exclude 
experimentally. It is not easy to detect y decay vertex, then we should take 2° decay effects in theoretical estimates. 
We find that we can simulate the decay effects by multiplying a factor 0.39 to C(q) - 1, if the pre-decay correlation in 
the EA or EE channel is small. For more serious studies, we need to take account of the EA and EE interaction and 
correlation in a given model of BB interaction consistently. 
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In the right panel of Fig. |4] we compare AA correlation using AA interactions under consideration. We have 
included the S° decay effects by the above mentioned simple prescription. We have not made the x 1 analysis, but 
the AA correlation data seems to favor fss2, NF(/? e = 0.50 fm) and some versions of NSC97 interactions. ND(R C — 
0.56 fm) and NSC89(m cllt = 1020 MeV) may be also allowed. These favored interactions are shown in the (1/ao, r e g) 
plane marked with open circles in Fig. [2] We conclude that AA interactions with 1/o.q < -0.8 fm -1 and r e ff > 3 fm 
are favored by the recent AA correlation preliminary data in high-energy heavy-ion collisions at RHIC by STAR 
collaboratio n 02811. These results are consistent with the analysis of the Nagara event, which is based on NSC97 
interactions lll9[l20ll . A more recent Nijmegen interaction, ESC08 J29ll . has a similar scattering parameters to fss2, 
and it is also in the above mentioned range. 
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Figure 5. H signal in AA correlation. The yield of H is assumed to be 20 times larger than the statistical model result. 



There are some more effects to be discussed for a more quantitative discussion on AA interactions. Another 
one is the coupling effects, AA «-» EN. Here we have included the coupling effects in a simple coupling potential, 
Vaa-s/vC'') = v cpi exp(-r 2 //? 2 ) with b — 1 fm. When the coupling is strong, the AA component is suppressed in the 
inner part of the relative wave functions. As a result, the coupling potential acts as repulsive interaction. We find that 
the coupling effects to EN channel are mild, as long as the coupling potential is not very strong, v cp i < 50 MeV. 

Finally, we discuss the existence of the H particle. We assume that the mass of H is 17 MeV above the AA 
threshold, and its width is 1.5 MeV. For the yield, the statistical model result is adopted as a reference value. Figure |5] 
shows the H signal which would be observed in the AA correlation, but with 20 times larger yield than the statistical 
model result. At present, we do not see any signal of H, but we need much more precision to conclude the existence 
of H from AA correlation unfortunately. This is because the number of background AA pair is large, then the signal- 
to-noise ratio is small for the H above the threshold. 



4. Summary 

We have discussed the exotic hadron size and hadron-hadron interactions in terms of the hadron yield and hadron- 
hadron correlation in heavy ion collisions. 

In the first part, we have demonstrated that heavy ion collisions may play a role of the hadron size ruler: In the 
framework of the coalescence model, a hadron with a large size would be produced more abundantly compared with 
the statistical model result. The mechanism of this enhancement is argued in a simple two body coalescence case, and 
we have found that coalescence favors hadrons whose shape in the phase space is similar to that of the source. If this 
coalescence mechanism also applies to other reactions, the yield of compact multi-quark states may be larger than 
the statistical model result in e + e~ reaction. Thus it will be valuable to investigate exotic hadron production in <? + e~ 
collisions. We have assumed that the source size is large enough in heavy-ion collisions, but it is necessary to take 
account of the finite size effects in e + e~, where the source size would be compatible with the hadron size. 

In the second part, we have discussed the AA interaction and its effects on the AA correlation, which is recently 
measured at RHIC f28ll . Based on a Gaussian source assumption and by considering the decay effects of 2° — > Ay, we 
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have compared the data with calculated results using several types of AA interaction. We find that the RHIC-STAR 
data favor the AA scattering parameters in the range 1 /oq < -0.8 frrT 1 and r e ff > 3 fm. These are consistent with 
the A A interaction parameters which reproduce the A A bond energy in AA He and recent A A interactions. The A A 
correlation data at low relative momenta seem to unfavor the existence of the bound H state containing significant 
AA component. In order to identify /rule out the resonance H state by using AA correlation data at higher relative 
momenta, we need more statistics because of the large continuum AA pair yield. For more serious estimate of the AA 
interaction, it is necessary to combine the feeddown, flow, and couple channel effects simultaneously. 
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